MG  2 2 1 978 


AEDC-TR-77-119 

ckj.Z 


INFRARED  SOURCE  TO  SIMULATE 
EARTH-LIMB  BACHRROUND 


Willis  L.  Clark 

ARO,  Inc.,  a Sverdrup  Corporation  Company 


VON  KARMAN  GAS  DYNAMICS  FACILITY 
ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
AIR  FORCE  SYSTEMS  COMMAND 
ARNOLD  AIR  FORCE  STATION,  TENNESSEE  37389 


August  1978 


Final  Report  for  Period  September  5 — 16,  1977 


Approved  for  public  release;  distribution  unlimited. 


Property  of  U.  S.  Air  Pore® 
AEOC  LIBRARY 
r 40600*77-0-0003 

Prepared  for 

ARNOLD  ENGINEERING  DEVELOPMENT  CENTER/DOTR 
ARNOLD  AIR  FORCE  STATION,  TENNESSEE  37389 


NOTICES 


When  U.  S.  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other 
than  a definitely  related  Government  procurement  operation,  the  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever,  and  the  fact  that  the  Government  may  have 
formulated,  furnished,  or  in  any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is 
not  to  be  regarded  by  implication  or  otherwise,  or  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell 
any  patented  invention  that  may  in  any  way  be  related  thereto. 


Qualified  users  may  obtain  copies  of  this  report  from  the  Defense  Documentation  Center. 


References  to  named  commerical  products  in  this  report  are  not  to  be  considered  in  any  sense 
as  an  indorsement  of  the  product  by  the  United  States  Air  Force  or  the  Government. 


This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is  releasable  to  the  National 
Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be  available  to  the  general  public, 
including  foreign  nations. 


APPROVAL  STATEMENT 

This  report  has  been  reviewed  and  approved. 

c . 

HERMAN  E.  SCOTT 

Project  Manager,  Research  Division 

Directorate  of  Test  Engineering 


Approved  for  publication: 
FOR  THE  COMMANDER 


ROBERT  W.  CROSSLEY,  Lt  Colonel,  USAF 
Acting  Director  of  Test  Engineering 
Deputy  for  Operations 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1 . REPORT  NUMBER 

AEDC-TR-77-119 

2.  GOVT  ACCESSION  NO. 

1 

3.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle)  5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 


INFRARED  SOURCE  TO  SIMULATE  Final  Report-September 

EARTH-LIMB  BACKGROUND  5 - 16,  1977 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR?*; 

Willis  L.  Clark,  ARO,  Inc. 

9 PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Arnold  Engineering  Development  Center/DOT 
Air  Force  Systems  Command 
Arnold  Air  Force  Station,  Tennessee  37389 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Arnold  Engineering  Development  Center/DOS 
Arnold  Air  Force  Station,  Tennessee  37389 

14.  MONITORING  AGENCY  NAME  & A D D R ESS(il  difleren  t from  Controlling  Office) 


15a.  DECLASSIFICATION7  DOWNGRADING 

schedule  n/a 

16.  DISTRIBUTION  STATEMENT  (of  this  Report)  " ~ ~ ™ “ — ” — ‘ 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  enfered  In  Block  20,  if  different  from  Report) 


18.  supplementary  NOTES 

Available  in  DDC 


19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

test  methods  altitude  space  environments 

earth  limb  surveillance 

simulator  gradients 

temperature  heat 

20.  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  identity  by  block  number) 

This  report  describes  the  preliminary  design  and  check  out 
of  an  earth-limb  simulator.  The  actual  simulator  will  operate 
in  a temperature  range  from  40°K  to  100°K  for  an  altitude  range 
of  60  to  115  km.  In  order  to  check  the  design,  temperatures 
were  uniformly  elevated  for  the  preliminary  testing. 


8-  CONTRACT  OR  GRANT  NUMBERfsJ 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  & WORK  UNIT  NUMBERS 

Program  Element  65807F 


12.  REPORT  DATE 

August  1978 

13.  NUMBER  OF  PAGES 

35 

15.  SECURITY  CLASS.*  (of  Ih Is  report) 


UNCLASSIFIED 


DD  , j an  ^73  1473 


EDITION  OF  1 NOV  65  IS  OBSOLETE 


UNCLASSIFIED 


AEDC-TR-77-1  19 


PREFACE 


The  work  reported  herein  was  conducted  by  the  Arnold 
Engineering  Development  Center  (AEDC) , Air  Force  Systems 
Command  (AFSC) , under  Program  Element  65807F.  The  results 
of  the  test  were  obtained  by  ARO,  Inc.,  AEDC  Division  (a 
Sverdrup  Corporation  Company) , operating  contractor  for  the 
AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  ARO 
Project  No.  V35S-B5A.  Herman  E.  Scott  was  the  Air  Force 
project  manager.  The  data  analysis  was  completed  on 
September  16,  1977,  and  the  manuscript  was  submitted  for 
publication  on  November  7,  1977. 
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1.0  INTRODUCTION 


A satellite-borne  surveillance  sensor  scans  from  outer 
space  to  the  hard  earth  in  its  search  for  a target.  In  such 
a scan,  the  sensor  encounters  a thermal  gradient  which  in- 
creases as  the  scan  approaches  earth.  This  thermal  gradient 
(which  starts  with  the  cold,  rarefied  gases  of  the  upper  atmos- 
phere and  continues  to  the  hot,  dense  gases  at  sea  level)  is 
referred  to  as  the  earth-limb  thermal  gradient.  The  earth- 
limb  gradient  and  a typical  satellite-borne  sensor  scan 
sequence  are  illustrated  in  Fig.  1 . 

The  calibration  of  such  a sensor  system  requires  the 
investigation  of  phenomena  directly  attributable  to  the  mission 
requirements  of  scanning  into  the  earth-limb  region.  A sensor 
viewing  a missile  through  the  earth-limb  region  requires  reason- 
able off-axis  rejection  to  separate  the  missile  from  the  back- 
ground radiation  of  the  earth  limb.  Should  the  mission  require 
the  sensor  to  scan  further  down  toward  hard  earth,  the  sensor  will 
go  into  full  saturation.  It  is  necessary  to  know  where  satura- 
tion will  occur  with  respect  to  the  viewing  altitude. 

If  the  sensor  (while  in  saturation)  starts  a scan  away 
from  hard  earth  toward  outer  space  it  will  remain  in  saturation 
for  some  finite  period  of  time.  The  recovery  time  is  of  utmost 
importance  for  a surveillance  satellite. 

An  earth-limb  simulator  concept  for  use  in  sensor  testing 
is  shown  in  Fig.  2.  The  thermal  gradient  of  the  earth-limb 
simulator  is  projected  by  the  imaging  mirror  onto  a diffuse 
screen  then  projected  as  a collimated  beam  to  fill  the  entrance 
aperture  of  the  sensor . 

The  Optical  Signatures  Code  (OSC)  is  a large  collection 
of  computer  programs  which  predict  the  radiation  emitted  by 
typical  targets  and  backgrounds  viewed  by  long  wavelength  infra- 
red (LWIR)  sensors  in  ballistic  missile  defense  missions.  The 
programs  in  the  OSC  were  developed  under  the  sponsorship  of  the 
Army's  Ballistic  Missile  Defense  Advanced  Technology  Center, 
and  some  are  now  operational  on  the  Arnold  Engineering  Develop- 
ment Center  (AEDC)  central  computer.  A plot  of  the  spectral 
radiance  at  100-km  altitude  obtained  from  the  OSC  is  presented 
in  Fig.  3.  The  spectral  radiance  data  from  the  Earth  Limb 
Signatures  Code  was  used  to  establish  the  radiation  requirements 
for  the  earth-limb  background  simulator. 
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2.0  DESIGN  AND  TEST  PROCEDURES 


A typical  LWIR  sensor  is  sensitive  to  radiation  in  the 
9-  to  11-ym  band  and  has  a 2.5-deg  field  of  view.  Installed 
in  the  AEDC  von  Karman  Gas  Dynamics  Facility  (VKF)  Aerospace 
Chamber  (7V)  (Chamber  7V) , and  looking  directly  at  the  23-in. 
collimating  mirror,  a 10-in.  by  10-in.  plate  is  required  at 
the  focal  point  of  the  mirror  to  fill  the  field  of  view  of 
the  sensor. 

The  radiance  values  integrated  over  the  9-  to  11-ym 
band  covering  the  appropriate  altitude  range  was  obtained 
from  the  Earth  Limb  Signatures  Code.  These  values  were  then 
substituted  into  Planck's  equation  to  calculate  the  black- 
body  temperature  required  to  produce  the  desired  radiance 
level  at  each  altitude.  The  temperatures  calculated  by  this 
procedure  were  in  the  30°K  to  100°K  range. 

This  study  resulted  in  the  design  of  a heated  plate  with 
a temperature  profile  to  match  the  temperatures  calculated 
above.  The  plate  was  tested  at  ambient  conditions  to  allow 
standard  temperature  measuring  techniques  and  to  check  calcu- 
lations used  in  the  mathematical  model.  The  calculated  tem- 
peratures were  elevated  to  the  32°C  to  60 °C  temperature  range, 
and  the  resulting  temperature  profile  is  shown  in  Fig.  4. 

Using  the  minimum  plate  dimensions  and  the  temperature 
profile,  a stainless  steel  plate  of  varying  thickness  which 
would  reproduce  the  required  temperature  profile  was  designed. 
The  design  required  the  derivation  of  the  appropriate  heat- 
transfer  equations  and  boundary  conditions  for  the  mathematical 
model.  The  mathematical  model  using  the  Runge-Kutta  integra- 
tion method  is  shown  in  Appendix  A. 

The  design  resulted  in  a plate  heated  on  both  ends  and 
varying  in  thickness  as  shown  in  Fig.  6.  The  calculated  heat 
input  was  10.9  watts  for  the  warm  end  and  5.03  watts  for  the 
cool  end.  These  inputs  assumed  the  plate  to  be  under  vacuum 
and  radiating  to  a 21 °C  heat  sink.  Figure  5 is  the  tempera- 
ture profile  of  the  plate  with  the  given  inputs  calculated 
from  the  method  shown  in  Appendix  A. 


3.0  RESULTS 


When  the  fabrication  of  the  plate  was  completed,  an  IR 
thermal  imaging  camera  was  used  to  acquire  the  temperature 
profiles  on  the  plate.  The  thermal  imaging  camera  is  a 
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television-like  scanning  system  which  displays  thermal  images 
on  a monitor.  The  camera  utilizes  a liquid-nitrogen-cooled 
Indium  Antimonide  detector  which  has  spectral  sensitivity  over 
the  2-  to  5.6-ym  range.  The  minimum  detectable  temperature 
difference  is  better  than  0.2°C  at  an  object  temperature  of 
30°C.  Frames  are  produced  at  a rate  of  16  per  second  which 
corresponds  to  a line  frequency  of  1,600  lines  per  second. 

The  images  produced  by  this  system  (by  remote  sensing  of  in- 
frared radiation)  provide  instant  thermal  profiles. 

The  data  acquisition  technique  was  by  adjustment  of  the 
isothermal  contour  lines  on  the  output  monitor.  The  controls 
provide  two  isothermal  contour  lines,  the  temperature  differ- 
ence being  measured  as  the  difference  between  the  two  iso- 
thermal contours.  In  Fig.  14,  the  two  vertical  lines  within 
the  green  color  band  are  the  two  isotherm  markers.  The  marker 
on  the  left  corresponds  to  the  isotherm  at  the  center  of  the 
plate  indicated  by  the  center  arrow.  The  marker  on  the  right 
corresponds  to  the  isothermal  contour  lines  indicated  by  the 
two  arrows  on  each  side  of  center.  Figures  14  through  13  are 
photographs  of  the  monitor  at  various  temperatures  and  are 
typical  for  the  full  range  of  data.  Figure  13  is  a photograph 
of  the  TV  monitor  output  without  the  isothermal  contour  lines. 

Figures  7,  8,  and  9 are  the  plots  of  the  temperature  data 
for  the  middle,  left,  and  right  sides  of  the  plate  along  with 
a plot  of  the  design  goal.  It  can  be  seen  from  these  figures 
that  the  general  profile  of  the  temperature  was  reproduced  by 
the  plate.  No  attempt  was  made  to  match  the  absolute  level 
between  the  required  and  acquired  data  because  of  the  long 
stabilization  time  inherent  in  the  plate.  More  input  heat  was 
required  than  was  calculated;  however,  the  test  was  conducted 
under  ambient  conditions  rather  than  under  vacuum  conditions 
as  used  for  design  calculations. 

Figures  10,  11,  and  12  are  various  views  of  the  plate. 

The  shape  of  the  plate  can  be  seen  in  Fig.  12  and  compared  to 
the  shape  shown  graphically  in  Fig.  6. 


4.0  CONCLUSIONS 


The  general  profile  of  the  earth  limb  was  reproduced  in 
a flat  plate  which  now  requires  only  the  modification  of  the 
input  power  to  be  used  as  a simulator  for  the  low  level  earth- 
limb  radiation.  Further  tests  should  be  conducted  to  deter- 
mine the  effect  of  various  background  temperatures,  the  effect 
of  vacuum  conditions  on  power  input  and  temperature  distribu- 
tion, and  possibly  the  effect  of  various  emissive  coatings. 
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From  the  data  obtained  in  this  check  out,  the  modeling 
calculations  were  adequate.  Some  further  refinement  may  be 
required  to  predict  the  input  energy  requirements;  however, 
it  will  require  study  in  a vacuum  environment  to  make  these 
final  adjustments  to  the  equations.  As  shown  in  the  data 
plots,  the  shape  of  the  required  temperature  profile  was 
simulated  by  the  plate.  Further  study  is  required  to  provide 
a plate  with  the  same  contour  but  with  temperatures  in  the 
30°K  to  100°K  range  as  required  for  sensor  calibration. 
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Figure  1.  LWIR  sensor  environment. 
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Figure  4.  Temperature-altitude  curve  for  earth  limb. 
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Figure  5.  Required  plate  temperature-plate  dimension. 
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Figure  6.  Calculated  plate  temperature-plate  dimension. 
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Figure  7.  Section  of  test  data 
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Figure  8.  Section  of  test  data  at  left  side  of  plate. 
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Figure  9.  Section  of  test  data  at  right  side  of  plate. 
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Front  view  of  simulator 
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Figure  12.  Side  view  of  simulator. 


Figure  11.  Rear  view  of  simulator. 
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Figure  13.  Thermal  image  of  simulator  without 
isothermal  contour  lines. 
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Figure  16.  Thermal  image  of  simulator  with  isothermal 
contours  set  at  = 10  C. 
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Figure  17.  Thermal  image  of  simulator  with  isothermal 
contours  set  at  AT  = 20  C. 
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APPENDIX  A 

EARTH-LIMB  HEAT  TRANSFER  FOR  CONTROLLED  TEMPERATURE 
DISTRIBUTION  OF  A RADIATING  SURFACE  IN  VACUUM 


For  a plate  radiating  on  two  sides  and  with  the  edges 
insulated,  the  heat  equation  can  include  variable  thickness, 
t,  as  follows,  assuming  one  side  is  flat. 


where  0 = slope  of  the  plate  contour  side 
= tan  dt/dx 

Using  the  Runge-Kutta  integration  technique  and  initial 
conditions  x 20,  T = Tq,  dT/dx  =0,  t = t , and  dt/dx  = dtQ/dx. 

A program  has  been  established  using  the  HP  9872  computer 
for  temperature  along  the  plate  with  any  given  T , T , t,  and 
dt/dx.  ° s 

Where  temperature  is  given  along  the  plate  and  thick- 
ness of  plate  required,  a Simpson's  rule  integration  of  total 
heat  can  be  used  as  follows 


Assume  small  increments  of  plate,  Ax,  and  variable  thickness  t. 


tl  = 


— 2 4 

a e Ax  Ts 

2 K AT 


" 1 <!  + coi-e> 
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6 = tan 


, t - t, 

-1  O 1 

Ax 


Then 


— 2 

n a e Ax 

+-  = y p 

n J K AT  s 


4 


'T  + T , \ 4 
n n+1 

V 2 T / 


- 1 


1 + sec  tan 


. (t  - t , \ 
-1  | n n-1 

V Ax  1 


—2 

CT  £ AX 

i 

2 K AT  Ts 


^n  + n-1 ) 4 

i 2T_  / 


- 1 

J 


1 + sec  tan 


, t - t J 
-In  n-1 


Ax 


Units:  Ax  = in. 


a = Btu/ftz/°R  K = Btu/ft 


2 4 

n a e Ax  T 

t = Z S 


'n  12  K(T  -T  , ) 

1 n n-1' 


- i 


1 - sec  tan 


a e Ax2  T 4 
s 

12  K(T  -T  . ) 
n n-17 


'T  + T , \ 4 
n n-1  ] 

, 2 T 


- 1 


1 - sec  tan 


-1 


ft  - t ^ 
, n n-11 

i AX  I 


Using  curve  fit  equations  for  temperature  and  assuming  parabolic 
form  for  points  near  zero  slope  to  expedite  and  simplify  the 
curve  fit. 


and 


T’ 


T 


AT+0 


dT 

dx 


4 
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where 


t = material  thickness  at  T 
n n 


and 


6 = tan  ■*"  dt/dx  = tan  ^ 


n 


- t 


n-1 


Ax 


9 is  assumed  to  be  zero  for  starting  the  integration  scheme 
and  hold  zero  until  t/x  was  significant  to  avoid  the  dilemma 
at  small  values  of  T'.  Also,  the  parabolic  form  of  curve  fit 
for  small  values  of  T'  was  matched  in  T and  T*  at  x = 1 given 
by  best  curve  fit.  Then  from  the  parabola 


x=l 

K 


KX 

T 


n 


x=l 


and 


T'  = nKX' 


n = T'/T 


n-1 


x=l 


(T'/T 

x=,  " |x=l 
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